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where  ^  is  the  surface  traction,  _u.  is  the  displacement, 

o 

is  the  portion  of  the  crack  generated  in  the  time  interval  [t  ,t], 
and  -£(t)  =  z  -  z_  with  z^  at  the  position  of  the  crack  tip  at 

t  _r  t  j- 1 

time  t.  o 

To  simplify  our  analysis,  we  avoid  geometrical  and  notational 
complications  by  limiting  our  discussion  to  edge  cracks  in  two- 
dimensional  bodies.  Also,  our  analysis  is  based  on  classical 
smoothness  hypotheses  and  therefore  in  applying  our  results  care 
must  be  taken  to  insure  that  the  underling  neighborhood  of  the 
crack  tip  is  free  of  shock  waves,  etc. 
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1.  Introduction  / 

it  is  the  purpose  of  this  paper  to^fve  a  unified  treatment 
of  the  dynamic  energy  release  rate,  (  t,J  for  a  sharp,  straight 
crack  in  a  hyperelastic  body  undergoing  finite  strain.  As  our 

■  '  Pf  J 
\  - 

main  result  we  decompose  'yt  into  the  usual  quasi-static  energy 
release  rate  plus  a  nonpositive  dynamic  contribution;  thus  for 
a  dynamic  solution  the  energy  release  rate  computed  using  the 
classical  quasi-static  formula  is  larger  than  the  actual  dynamic 
energy  release  rate.V^We  also  present  what  are  apparently  the 
first  proofs  (within  the  dynamic  theory)  of  the  well  known  relations 


e<to) 


V 


tot 


where  s  is  the  surface  traction,  u  is  the  displacement, 

t  is  the  portion  of  the  crack  generated  in  the  time  interval 
o 

(t  ,t] ,  and  i(t)  =  z  -  z  with  z  the  position  of  the  crack 
tip  at  time  t. 
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To  simplify  our  analysis,  we  avoid  geometrical  and 
notational  complications  by  limiting  our  discussion  to  edge 
cracks  in  two-dimensional  bodies.  Also,  our  analysis  is  based 
on  classical  smoothness  hypotheses^-  and  therefore  in  applying 
our  results  care  must  be  taken  to  insure  that  the  underlying 
neighborhood  of  the  crack  tip  is  free  of  shock  waves,  etc. 


Notation.  Light-face  letters  indicate  scalars;  bold- face 
lower  case  letters  indicate  vectors  (in  JR  )  ;  bold-face  upper 

case  letters  indicate  second-order  tensors  (linear  transformations 

2  2  T  2 

from  JR  into  3*  )y  £  is  the  transpose  of  A;  A»B  = 

div  s  is  the  vector  with  components  dS^/dx..?  ?u  is  the  tensor 

with  components  Buj/Bx^ ;  V7u  is  the  third-order  tensor  with 

components  3  u^/Bx^Bx^;  a  superposed  dot  denotes  differentiation 

with  respect  to  time;  Lp(ft)  is  the  class  of  all  functions  <p 

on  ft  with  |tp)p  integrable  on  ft. 


^Such  hypotheses  are  tacit  in  most  other  studies  of  this  type. 

An  exception  is  Freund  [1977] . 

Here  we  use  standard  indicial  notation  and  cartesian  coordinates. 
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2.  Basic  equations . 

To  fix  notation  we  consider  first  a  two-dimensional  regular 

2 

body  B,  which  we  identify  with  the  regular  region  of  ]R 
it  occupies  in  a  fixed  reference  configuration.  We  assume  that 
the  body  is  hyperelastic, ^  so  that  the  displacement  u(x,t) , 
the  (Piola-Kirchhoff )  stress  S(x,t),  and  the  stored  energy 
w(x,t)  obey  the  energy  equation 

w  =  £-Su  (2.1) 

and  the  equation  of  motion 

div  S  =  pu  (2.2) 

with  p  >  0  the  density  in  the  reference  configuration.  We 
assume  throughout  that  p  is  constant . 

The  above  equations  are  appropriate  to  both  the  finite 
and  infinitesimal  theories  of  elasticity,  in  the  infinitesimal 
theory  S  is  symmetric  and  w  quadratic,  but  these  restrictions 
are  not  relevant  to  most  of  what  follows. 


H»ith  the  exception  of  Section  6,  our  analysis  is  valid  for  more 
general  materials  provided  one  uses  (2.1)  as  the  definition  of 


w. 


We  consider  a  fixed  time 


3  .  Mathematical  preliminaries . 


interval  [0,T]  .  We  assume  that1 

B  contains  a  straight  edge  crack 

Ct«  The  tip  of  the  crack  at  time 

t  is  denoted  by  z^;  we  assume 
2 

that  z.  is  C  in  t  with 
~t 

velocity 

£(t)  =  it  5t  *  2  t3-1) 


and  that 

o 

z,  eB 
~t 

o 

for  0  <£  t  <£  T,  where  B  is  the  interior  of  B. 

The  fields  cp(x,t)  of  interest  will  be  defined  at  each  x 
in  B\Ct  and  each  te [0,T] .  A  field  of  this  type  is  a  Cn 
fracture  field  (n  J>  0  an  integer)  if: 

(i)  the  derivatives  of  tp  of  order  <£  n  exist  away 
from  the  crack? 

(ii)  «p  and  its  derivatives  of  order  <£  n  are  continuous 

away  from  the  crack  and,  except  at  the  tip,  are  continuous 
up  to  the  crack  from  either  side. 


We  write 

(B ) 


It  is  important  to  note  that  B  here  need  not  be  the  entire 
body,  but  rather  an  arbitrarily  small  neighborhood  of  the  tip 
(cf .  the  remark  preceding  Theorem  1) . 


c 


if  <p(*,t)eLP(8)  at  each  te[0,T].  If  <peLP(B) ,  then  given 
any  one-parameter  family  8^  (5  >  0)  of  regular  subregions 

of  B  with  area(8\Bg)  -»  0  as  6  -*  0, 

J  |<p(x,t)  |p  da  -*  J  I(p(x,t)  |p  da 

8.  B 

0 

as  6  -*  0;  when  this  limit  is  uniform  in  t€[0,T],  we  say 
that  <peLp(B)  uniformly  in  time .  Analogous  interpretations 
apply  to  the  assertions  «peLp(SB),  tpeLP(SB)  uniformly  in  time,  etc. 

Consider  now  <p(x,t)  as  a  function  cp(zt+r,t)  of  t  and 
the  position  vector  r  from  the  tip.  We  let  denote  the 
derivative  of  this  function  with  respect  to  t  holding  £ 
fixed;  thus,  by  (3.1), 

{ p^=  <p  +  V  (p  (3.2) 

c 

with 

?C<P  =  vjp's 

i 

j 

the  directional  derivative  of  o  in  the  direction  c.  For  a 
vector  field  u,  u^  is  defined  in  the  same  manner,  except  that 
now 

V*  = 

Since  Ct  c:  CT  for  0  £  t  <£  T,  each  Cn  fracture  field 
is  of  class  Cn  on  the  cartesian  product  (B\CT)  X  [0,T]  .  The 
next  two  lemmas  give  certain  important  identities  for  functions 
of  this  type. 


L 


5a 


Lemma 


1.  Let  <p  be  a  C"*"  scalar  fie  Id  on~*~  (B\C^,)  x  [0,T]  . 


Assume  that  <p,<peL  (B)  uniformly  in  time.  Then 


J  cp  da 


B 

is  a  function  of  time  and 


&  j  ,  da  =  J'  4  da. 
B  B 


LOf  course,  <p  may  have  singularities  on  CT , 


(3.3) 


(3.4) 
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Proof . 

Let  ft. 

0 

_ mm _ 

L- 

I 

St 

V 

J4  -  B\»4’ 

so  that  < 

of  a  portion  of  SB,  two  lilies 
parallel  to  C^,  and  a  straight 
end  face  of  length  6  perpendicular 
to  these  lines  and  such  that  C„ 


is  contained  in  ft^ 
on  Bj  x  [0,T]  .  Then 

f* 


Let 


^  I  <P  da  =  j  <P  da- 
B.  B. 


(3.5) 


Moreover,  as  tpeL  (B)  uniformly  in  time. 


Jtpda-*J<pcla  (3.6) 

Bft  B 

as  6  -*  0,  uniformly  in  time.  Thus,  since 

J  <p  da  -*  J  <p  da,  (3.7) 

*4  0 

(3.4)  holds.  Further,  since  the  left  side  of  (3.6)  is  continuous 
and  the  convergence  uniform,  the  right  side  must  also  be  continuous; 
hence  (3.3)  is  C1  in  time.  □ 


Henceforth,  in  boundary  integrals  the  letter  n  will 


always  designate  the  outward  unit  normal. 


n 


Lemma  2 .  Let  <p,  and  tu  be  scalar  fields  on 

(fi\CT)  X  [0,T] .  Assume  that: 

(i)  <p  and  u)  are  C*  and 


«p  =  il)  +  ?c«>; 

(ii)  q>,jJ)eL^(8)  and  uoeL^{98),  all  uniformly  in  time. 


Then 


J  cp  da 

8 


is  in  time  and 


•“  f{pda=JjJ)da+J  u)  c-n  dAf. 
8  8  50 


Proof .  Let  and  8fi  be  as  in  the  previous  proof.  On 

the  upper  and  lower  horizontal  portions  of  98  c*n  =  0,  where 
n  is  the  outward  unit  normal  to  98^.  Also,  by  the  continuity  of 
in  away  from  the  crack,  the  integral  of  «ic*n  over  the  vertical 
right  end  face  of  Rg  tends  to  zero  as  6  -*  0  uniformly  in  time. 
Thus  (ii)  and  the  divergence  theorem  imply  that 


ujc  *n  cU>  = 
J  +*  ** 

98 


lim 

t-*0 


I 


ujc.n 


98  6 


dd/  = 


lim  7  it)  da; 

6-*°  J  C 
B6 


hence  (i)  and  (ii)  yield 


f  ii)  da  +  f  luc-n  dAS= 
J  j  **»  ~ 

8  98 


lim  («i)+V  uu)da 

6-*°  J  C 
8  t 


lim 

6-»0 


J  * 


8, 


da , 


and  this  limit  is  uniform  in  time.  But  (3.5)  and  (3.7)  hold  in 
the  present  circumstances .  Hence 


J  tp  da  =  lim  J  9  da  =  j  ip  da  +  J  uic*n  d A/, 

B  ^Bg  B  aft 

and,  as  in  the  proof  of  Lemma  1,  the  left  side  is  continuous  in 

time .  1  I 

-s.  In  the  next  lemma,  and  in 

/  N.  the  sequel,  fig  =  Clg(t)  is  the 

f  P  J  J  disc  of  radius  6  centered  at 

\  J  the  crack  tip,  z 

Lemma  3.  Let  f  be  a.  C^"  fracture  field  with  div  f€L^(ft)  . 
Assume  further  that 


lim  f  f  *n  6m  -  0, 

6-0  J,  ~  ~ 


-  0  as 


and  that,  for  e  a  unit  vector  normal  to  C.  ,  f(x.t)  *e 

■■■■"  '  -  ' -  '  '  .  1  ■■  —  ■£_  ?  ***'*s*7  ' 

x  approaches  Ct\{zt)  from  either  s ide .  Then 

f*n  cW  =  div  f  da. 

J  ~  ~  j  ~ 

bb  ft 


Proof .  Simply  apply  the  divergence  theorem  to  the  region 
B\n.  and  let  6—0.  I  1 
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4.  Basic  assumptions .  Energy  release  rate . 

We  begin  by  stating  our  assumptions  concerning  the  fields 
u(x»t)»  SLOSjt)*  and  w(x,t)  .  Here  and  in  what  follows  the 
field  s  in  a  boundary  integral  will  always  denote  the  surface 
traction 


s  =  Sn. 


(4.1) 


•  3  1 

(A.)  u  is  a  C  fracture  field;  s  and  w  are  C 

J.  ~  »» 

fracture  fields;  u,  S,  and  w  obey  (2.1)  and  (2.2)  away  from 
the  crack. 


(A  )  For  f  =  u,  u#,  u1,  u'"j  f  is  bounded;  S,VfeL2(B) 
uniformly  in  time;  WjW^VVfeL1 (B)  uniformly  in  time. 

(A3)  The  surface  traction  vanishes  on  the  crack;1  that  is, 
if  e  is  a  unit  vector  normal  to  C.  ,  then  S(x,t)e  -*  o  as  x 
approaches  C^Xfz^}  from  either  side. 

(A^)  Given  any  bounded  vector  field  v  on  B, 


lim  s-v  d JU-  0. 

«  _  ,J  *+  ** 

6-»0. 


'an, 


Assumption  (A^)  asserts  that,  for  bounded  "velocity  fields", 
the  virtual  power  vanishes  at  the  tip.  By  taking  v  equal  to  a 
rigid  velocity  field  one  immediately  concludes  that  (A^)  rules 
out  the  possibility  of  a  concentrated  force  or  moment  at  the  tip. 


Our  results  extend  trivially  to  the  case  in  which  tractions  are 
prescribed  over  CG,  the  initial  configuration  of  the  crack;  we 
simply  replace  integrals  over  aB  involving  s  by  corresponding 
integrals  over  aB  +  CQ.  Here  an  integral  over  CQ  has  the  obvious 
meaning  in  terms  of  integrals  over  the  "two  faces"  of  Co  (cf. 

(6.4))  . 

A 
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Finally,  (A^)  is  implied  by  the  somewhat  more  stringent 
assumption 


lim 

6-*0 


Is |d/d/  =  0. 


The  function  t  on  [0,T]  defined  by 

C=Js.u<U-^J  (w+k) da 
3ft  ft 

is  called  the  dynamic  energy  release  rate .  Here 


z  ~ 


is  the  kinetic  energy  per  unit  volume. 


(4.2) 


To  see  that  &  is  well  defined,  note  first  that,  by  (A2) 
and  the  identity 

u  =  u7-  v  u,  (4.3) 

~  ~  Cr' 

k  satisfies 

k,keL1(B)  uniformly  in  time. 

Thus  the  existence  of  the  second  term^  in  (4.2)  follows  from  the 
equation 

(w  +  k)  *  =  (w  +  k/  -  (w  +  k) 
and  Lemma  2  with  ep=-tu  =  w  +  k,  $  =  (w  +  k)^ 

■^The  first  term  is  well  defined. 
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Let  ft  be  a  regular  subregion  of  ft .  We  say  that  ft 
surrounds  the  tip  at  time  t  if  zteft  and  5ft  intersects 
Ct  at  only  one  point. 


Remark.  Let  ft  surround  the  tip  at  time  t_ .  Then  in 
-  o 

a  sufficiently  small  neighborhood  of  t  =  tQ,  (2.1),  (2.2), 
and  the  divergence  theorem  yield  the  energy  equation 

J  J  (w+k)da, 

3(B\ft)  B\R 

and  we  can  rewrite  (4.2)  in  the  form 


e  = 


aft 


*  Q  ' 

S-U  cU-  -  j  (w+k)  da  . 


ft 


Thus  6  is  intrinsic  to  the  crack  tip;  that  is,  the  definition 
of  t  _is  independent  of  the  region  ft  surrounding  the  tip. 


Our  first  result  gives  an  alternative  formula  for  C  in 
terms  of  a  boundary  integral  over  a  region  which  shrinks  to  the 
crack  tip,  z^.  With  this  in  mind  we  give  the  following  definition. 
Let  4  be  a  scalar-valued  set  function  defined  on  the  class  of 
all  regular  subregions  of  B.  Let  $QelR  .  We  write 

<f>  =  lim  $(ft) 

°  ft->z. 

'-t 

if  given  any  6  >  0  there  is  a  6  >  0  such  that 

\%  -  $(ft)  |  <  t 


i 
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for  every  region  ft  which  surrounds  the  tip  at  time  t  and 
has  area  less  than  5  . 


Theorem  1.  At  each  time  t, 


6 ( t)  =  lim 
ft-*z 


J  (<w 

t  3ft 


d  M 


(4 .4) 1 


Proof.  By  (2.2), 


div(gTu)  =  S-7u '+  pu.u' 


(4.5) 


and,  since 


u  =  u* -  27  u/  +  7ju  -  7. u,  (4.6) 

(v  c~  C"*  C~ 

(A2)  implies  that  the  right  side  of  (4.5)  belongs  to  L^fB)  .  Thus 
(A^) ,  (A^) ,  and  Lemma  3  imply  that 

r  s  •  u^  d-4'  =  I  (S  *7u^+  pu-u)da.  (4.7) 

3B  B 


Next,  by  (3.2),  (4.6),  and  (A2) , 


w  =  w  -  7  w, 
c  ’ 


k  =  K  - 


’c<2 


I’oHl2’ 


with 


K 


pu-^-  p^v!' -  27  u^  -  7^  -7  ue  L1(B). 


(4.8) 


Thus  Lemma  2  with 


<p  =  w  +  k, 


i5  =  v/+  K , 


U) 


2 


1  • 

For  a  linear  elastic  solid  a  similar  relation  (containing  -jj 
in  place  of  7^)  appears,  without  proof,  as  Eqt .  (13)  of  Atkinson 
and  Eshelby  [1968],  and,  with  a  sketch  of  a  proof,  as  Eqt.  (13) 
of  Freund  [1972] . 
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yields 

~  f  (w+k) da  =  J  (w^-K)da  -  J  (w  +  cU/, 

8  B  dB 

and  this  relation,  (4.2),  (4.3),  and  (4.7)  imply 

£  =  $(B)  -  |  (w^  +  K  -  S*7i/-  pu  *u)  da ,  (4.9) 

8 

where  £>(ft)  designates  the  integral  in  (4.4). 

Now  let  ft  surround  the  tip  at  time  t,  where  t  is 
fixed.  Then  the  remark  preceding  Theorem  1  implies  that  B 
in  (4.9)  can  be  replaced  by  ft: 

6  s  *(ft)  -  f  (w/+  K  -  5*7^-  pii-u^da. 
ft 

Since  the  above  integrand  belongs  to  L^(B),  (cf.  (A2)  and 
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(4.8)),  if  we  let  the  area  of  R  approach  zero,  we  arrive 
at  the  desired  result  (4.4)  .  1  I 


Remark .  It  is  important  to  note  that  (by  definition)  the  limit 
R  -*  z  in  (4.4)  need  not  be  confined  to  regions  R  whose 
diameters  tend  to  zero.  The  chief  requirement  is  that  the  area 
of  R  approach  zero.  Thus,  if  Rj  is  the  one-parameter  family 

of  regions  shown  in  the  figure, 
where  each  R^  is  a  rectangular 
region  which  surrounds  the  tip 
and  has  height  6,  width 
independent  of  6 ,  and  horizontal 
sides  parallel  to  the  crack,  then 

e(t)  =  lim  [  f  (w  +  -  ST^ou)  -n  d^; 

6_K)aR. 

0 

or,  since  £•#  vanishes  on  the  horizontal  portions  of  dR^ , 

while  the  integrals  over  the  vertical  portions  tend  to  zero  as 

6  0, 

£(t)  =  -  lim  I*  u)  -n  d^, 

6  -*0  J  ~  ~ 


=  -  lim  c»  [  VuTs  d 4J. 

6  ~*0  J 


Henceforth,  t  is  a  fixed  time  in  [0,T]  .  For  convenience, 
we  write 


(_c 1_\ 

\dt/t 


(4.10) 


for  the  derivative  with  respect  to  t  at  t  =  t  .  Also,  given 


a  fracture  field  <p,  we  write  for  the  function  on  B\Ct 

defined  by 


<Pt(x)  =  tp(x,t)  . 


Theorem  2 . 

t 

e<to>  =  felt  J  j  +  p<i4-3*>  •»*)«'  da.  (4.U) 

'  *  O  0  . 


Proof .  Let 

t 

3(t)  =  [  j*  s*u  d4/  dk  -  I  (w+k)da, 

t  3B  R 

o 

so  that,  by  (4.2), 

e  =  3. 

The  first  term  in  (4.12)  is  equal  to 
t  t 

j  J  *t’^dW  J  UA~St)  <4*14> 

t  3fl  t  30 

o  o 

^Cf .  Rice  [1965],  Eqt .  (5).  Rice  considers  infinitesimal  strains 
and  time- independent  applied  surface  tractions. 


(4.12) 


(4.13) 
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Since 


Aj  J  «Sx-£t>  -Sx  ^  ax  ■  -  J  J  it'iy,  <U/dx- 


t  5B 
o 


t  aB 

o 


(4.10)  applied  to  the  second  term  in  (4.14)  vanishes.  Thus 
t  t 

(lt)t  J  I  s-n  a*  ax  ■  (gj)t  J  j  St -ax  (^)t  J  ><U-, 


t  3B 
o 


dB  t 


and,  by  (A2)-(A4),  (2.2),  (4.1),  and  the  divergence  theorem 
(Lemma  3),  this  equals 

(dit)t  J  tet*(^f7St  >  +  Piit,(Ht-Ht  )}da;  (4*15) 

'  '  o  B  o  o 

hence  (4.13)  yields 

e(to}  =  fe)  t  J  tet^f^t  }  +  pBt'(iit^t  }  "  (w+M)da.(4.16) 

'  /  O  a  O  O 


Clear ly. 


>  “  J  £t,7i? 


for  points  x  not  on  .  Thus,  since  is  a  set  of  (area) 

measure  zero  in  B,  both  sides  of  this  equation  may  be  integrated 


over  B  to  give 


J  St*(7JSt-7^t  )da  =  J  J  £t’7^  dX  da* 
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Similar ly, 

t 

I  p£t*(Bt-JSt  )da  =  I  J  pSt-JSx  dx  da> 

8  °  B  t 

o 

and,  using  (2.1)  and  the  definition  of  k, 

t  t 

^  J  (OT-k)da  -  ^  J  J  (wk)dX  da  =  gj  J  J  +  pj^-u^dX  da. 

B  B  t  ft  t 

o  o 

The  last  three  results  and  (4.16)  yield  (4.11)  .  I  I 


Remark .  It  is  clear  from  the  remark  preceding  Theorem  1 
that  6  in  (4.11)  can  be  replaced  by  any  region  &  which 
surrounds  the  tip  at  time  tQ. 


There  is  much  confusion  in  the  literature  concerning  results  of 
this  type  (cf.  Gradin's  [1979]  correction  of  a  misconception  of 
Luxmoore  and  Morgan  [1977]  in  the  quasi-static  theory) . 
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5.  Partition  of  C. 

Theorem  2  allows  us  to  write 


£(t0)  =  U(tQ)  +  K(to) 


with 


(5.1) 


X  represents  the  contribution  of  inertial  effects  to  the  energy 
release  rate,  while  U  gives  the  value  the  energy  release  rate 
would  have  were  this  dynamic  contribution  neglected. 

For  a  linear  elastic  material,  with  positive  semi-definite 
elasticity  tensor,  or,  more  generally,  for  a  material  which  is 
stable  in  the  sense  of  Drucker  [1964] , 


U  2  0  • 

The  next  theorem,  which  is  our  main  result,  shows  that,  to  the 
contrary, 

X  £  0; 


hence  the  effect  of  inertia  is  to  reduce  the  energy  release  rate 
from  the  value  obtained  using  the  quasi-static  formula  (5.1)^. 


i 
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Theorem  3 .  Assume  that  u  is  continuous  at  t  in  the 

■  —  ii  —  —  -  -Tii  1  '  "  ' T  D  “  1  “  ■  ‘  ■ 11  1 


sense  that 


lim  sup  |u(x,t)  -  u(x,t) |  =  0 

t-»t  xefl  ° 

o  ~ 


(5.2) 


Then 


K(t0) 


*  J  ■§  <St-St  >da’ 


-  J  2  )-7o(t)(Bt-St  )da> 

"  o  o 


tjt  r  i  ivC(t)(Ht-Ht  >i  da* 


(5.3) 


hence  K  (t  )  <[  0 


Proof .  S ince 


I  <3t-St  >  dA  =  ^t-Bt  >  •  tet-2t  ]  > 


we  can  write  (5.1) 2  as  the  right  side  of  (5.3)^  plus  p/2  times 


(^)t  I  "i  da- 


(5.4) 


where 


tPl<x,t)  =  J  (ut  +  ut  -  »gA).S>  dA. 

‘o 

Thus,  to  establish  (5.3)^  it  suffices  to  show  that  (5.4)  vanishes 


For  convenience,  let 

3  =  Ht>  uo  =  ut  »  3  =  JJ-uo,  £0  =  c(tQ)  . 

A  simple  calculation  then  shows  that 

»  »*2  ®  2\ 

*1  =  <3  +  So5  *3  -  (JS  “  Jio)  • 

Next,  in  view  of  the  identities  (4.3),  (4.6),  and 


u'"  -  37  u" 
~  o~* 


we  have 


3?au/-  7«u  +  37^u/+  37  7 ai  -  73u, 
&•  c~  o~  C  O' 


(«PX-«P2)  =  ll)  +  7cu>, 


where 


»2  =  g.(^  -  ^)Uo, 

o 

ij)  =  (u'"  -  37  u"  +  37^u/-  37  «u^  +  37  7.u  -  7..u} 'Q  -  (li  _  )  .u^ 

~  c~  o~  c~  cc~  3  'Js  Ho'  H 

+  (u"  -  27  U/-  7,U  -  (u"  -  27  l/  -  7.  u  ) } *7  U 

O'  ~o  c  ~o  c  ~o  c~ 

o  o 

+  2q  •  7 . 7  u  +  l/-(7^  -  72  )u  , 

^  c  c^o  **  c  c  /  9 

o 

■  “  -s-^s  +  V’oa- 

It  is  easy  to  verify  that  ?p  =  tp^  -  <p2,  i!>,  and  iu  obey  the  hypotheses 

of  Lemma  2;  hence 

lit  J  ('Pi-tP2^da  =  J  ^  da  +  f  oic.n  dA4 

B  B  38 


¥ 
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and  since 


it  follows  that 


^(x,t  )  =  Ul(x,t  )  =  0, 


[dt)t  J  2)da  “  °* 


(5.5) 


Next, 


J  cp2da  =0  at  t  =*  t  ; 


hence 


where 


(^)to  J  *2da  = 

8  o 


*(tQ,t) 


t^r  I  <p2(£'t,dV 


Let  e  =  c/  |c  | .  Then 


q>,  =  o.  (c2-c2) V2u  , 
l  a  ~o  ~  e~o 


so  that 


2  2 
c  -c 


l^o't)  I  ^  (sup  |u-^i>RbHJ  i^0ida^ 


and,  toy  (5.2),  this  tends  to  zero  as  t  -♦  t  ,  since  V2u  eL1^) 

o  e~o  ' 


2  2 

c  -c  (c  -c) • (c  +c) 

<~o  w  '^*o  v~o  ~ ' 

t-t  "  t-t_ 


2  c  *c  . 
'~o  ~o 


Thus 


(5.6) 


and  (5.5)  implies  that  (5.4)  is  zero;  hence  (5.3)-^  holds. 
Next,  define 


*3  =  CHL>  *2» 


so  that  (5.3)^  becomes 


K(to»  •  -§fe)t  J  «3d“- 

°  B 


A  simple  calculation  shows  that 


cp3  =  v  +  <3.7^  -  <p2, 

y  -  3* (a*  -  2*02'-  "as  -  <a;  -  2%  s'a  -  »4  s0>)>  <5-7> 

o  o 

.  1 

and,  by  >  y,veL  (B)  uniformly  in  time.  Thus  y  obeys  the 

hypotheses  of  Lemma  1  and 


A  J  *  da  -  J  y  da, 

B  B 


But  v(x,tQ)  =  0  for  almost  every  xeB;  hence 


(at)t0  J 


■y  da  =  0, 


and,  by  (5.6)  and  (5.7)  j,,  (5.3)2  holds 


22 


and  (5. 3) ^  follows  from  (5.3)2* 

Finally,  since  the  integral  in  (5.3)3  vanishes  at  t  =  0 
and  is  ^  0  otherwise,  K (^q)  _<  0 .  I  1 

Remark .  It  should  be  emphasized  that  U  and  K,  separately, 
do  not  generally  represent  rates  of  change  of  strain  energy 
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w  = 

j*  w  da 

8 

and  kinetic  energy* 

K  = 

J  k  da. 

6 

Indeed,  it  is  clear  from  (5.1) 

that 

u(to)  =  -W<t0)  + 

l)t  J  )da> 

'  °  e 

(5.9) 

K(t0)  =  -K(tc)  f 

;)t  J  )da> 

°  8  ° 

When  the  power  supplied  by  the 

environment  vanishes,  i.e., 

when 

[  s-u 

J  w  w 

o 

II 

(5.10) 

38 

the  last  two  terms  in  (5.9)  sum  to  zero  (cf.  (4.15)),  but 

are 

generally  not  individually  zero,  so  that  while 

h  +  K 

=  -  W  -  K, 

in  general 

u  ¥  -w. 

X  ^  -K. 

This  is  in  contrast  to  the  quasi-static  theory  for  which  U  =  -W 
when  (5.10)  holds. 


*Cf.  the  theoretical  studies  of  Hahn,  Hoagland,  Rosenfield,  and 
Sejnoha  [1974],  Freund  [1977],  and  Popelar  and  Gehlen  [1979], 
where  actual  curves  are  given  for  strain  energy  and  kinetic 
energy  as  functions  of  crack  length. 
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Remark .  The  definitions  (5.1)  are  intrinsic  to  the  tip; 
that  is,  (5.1)  are  invariant  when  ft  is  replaced  by  an 
arbitrary  region  ft  surrounding  the  tip.  Indeed,  if  we 
replace  each  integral  over  ft  in  (5.1)  by  an  integral  over  ft 
plus  an  integral  over  ft\ft,  and  interchange  this  integral  over 
B\ft  with  the  time  derivative  at  t  =  tQ,  we  find  that  the  terms 
involving  B\ft  vanish . 

In  view  of  this  remark  and  the  remark  preceding  Theorem  1, 

C,  U,  and  K  are  possibly  more  closely  related  to  the  dynamical 

•  •  •  • 

behavior  of  the  crack  than  W  and  K,  since  W  and  K  are 
generally  not  intrinsic  to  the  tip. 


25 


6 .  The  energy  release  rate  for  a  linear  elastic  material . 

We  now  restrict  our  attention  to  the  linear  theory  for  which 

S  =  £7u,  w  =  ^  7U-C7U  (6.1) 

with  C,  the  elasticity  tensor .  a  symmetric  linear  transformation 
(at  each  point  of  B) .  Of  course,  symmetry  is  the  requirement 
that 

A  *CB  =  B -CA 


for  all  second-order  tensors  A  and  jg, 


Theorem  4. 


h(to) 


=  f  (S . -S.  )*(7u.-tfu  ) da, 

2  \dty  t  J  ~t  ~t  ~t  ~t 
'  '  °  B  o  o 

=  1^3 t-^t  )*£(^t-^t  )da'  (6*2) 

°  fl  ° 


G(t  )  =  4  ’  .  f  (  (S^-S.  )  •  (?u  -Vu.  )  +  p(ll.-u.  )  •  (u. -u.  ) 

'  o  2  Vdt,  t  J  1  '~t  ~t  ~t  ~t  —t  ~t  ~t 


B 


Proof .  As  in  the  first  few  steps  of  the  proof  of  Theorem  3, 
we  can  decompose  U (tQ)  as  follows: 


Mto)  =  ^ 


-  ife)t  I J  <St-St )  )da  +  J  J  (st+^t  ■ -2^>  dx  4 

°  B  °  °  B  to  °  J 


Since  C  is  symmetric, 


St.?ut  =  gt  *Vut, 
o  o 


XPor  (6.2^  cf.  Buechner  [1958],  Eq.  (15);  Rice  [1965],  Eq .  (19). 
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and  thus,  by  (2.1)  and  (6.1), 
t 

J  <£t+&t  dx  -  <£t+5t  >  -  o, 

£  O  o  o  o 

o 

which  yields  (6.2)^  j*  The  result  (6.2)^  is  a  direct  consequence 
of  (5 . 3)  1  and  (6.2)  1#  □ 


For  t  >  t  let 
o 


ct  t  =  ctxct 
o  o 

denote  the  portion  of  the  crack  between  and  z  . 

~t0  ~t 

statement  of  the  next  theorem  we  write 


In  the 


J  St  *2t  ^ 


(6.3) 


"tot 


for  the  integral  of  s  *u.  over  the  "two  faces"  of  0  .  ;  that  is, 

~T'  t_t 

.  .  °  O 

writing 

u^fxjt)  *  lim  u(x  +  6n+,t) , 

6-0+ 

si(x,t)  =  lira  S(x  +  6n+,t)r£, 

6-0+ 

for  all  xeCT  ft  zt)  ,  where  rp  =  +  e  with  e  a  unit  vector 

perpendicular  to  the  crack,  then  (6.3)  is  defined  to  be 


J  st  ^ d* +  J  s;  *#. 


(6.4) 


'V 


V= 


Also,  in  what  follows  the  notation  (4.10)  designates  the  right-hand 


derivative  at  t  . 

o 


27 


+  ,1 


Theorem  5.  Assume  that  sr  eL  (C„)  .  Then 

r  *  -=”  _  X 


=  -  ifdtk  I  5t  -JSt  w 

'  '  O  -  o 


V 


(6.5) 


lf_d) 
-  -  2  1  - 


(A)t  J 


'V 


where  t(t)  =  z,  -  z.  . 

o 


Proof .  By  (6.2)^, 


e(to> 


=  -4—)  r 

4dtM 


V  da. 


where 


(6.6) 


v  =  ^-Eo5  -*2  +  P(jj-)io)  -3 


(6.7) 


(Recall  that  g  =  u-uQ.)  Consider  the  region  &6  =  fi\{n6(to)  U  Og(t)) 


and  fix  te(tQ,T].  By  (2.2), 


y  =  div{  (S-S0)Tg)  . 

We  therefore  conclude  from  the 
divergence  theorem  and  (A.j) 
that 


For  the  quasi-static  theory  (6.5) 1  was  derived  by  Bueckner  [1958] 
(see  also  Rice  [1965]),  while  (6.5)2  is  essentially  contained  in  the 
work  of  Irwin  [1957,  1958] .  For  the  dynamic  theory  a  heuristic 
argument  in  support  of  (6.5)2  was  given  by  Erdogan  [1968].  The 
observation  that  (6.5)  require  separate  proof  in  the  dynamic  theory, 
and  that  such  a  proof  is  probably  non-trivial,  is  due  to  Freund 
[1972],  who  notes  that  (6.5)  yield  "the  correct  result  for  all 
problems  for  which  a  solution  is  known". 
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j  y  da  - 

B6 

J  (s-So>  '3  <3*'  ~  J  <2-S0>  '3  44'"  J  <£-So>-adi/+  J  <S-£o> 
9e  »0t(to)  an{(t)  B6net 

o 

where  the  last  integral  has  a  meaning  analogous  to  (6.4).  Clearly, 

J  s-a  da/-,  o 

a06 (‘o’ 

as  6  -♦  0,  because  both  S  and  q  are  bounded  near  z.  . 

~  ~  o 

Similar ly, 

J  ~o  *2  0  * 

an6(t)~° 


On  the  other  hand,  since  g  is  bounded,  (A^)  implies  that 

So *2  a*-*  °> 


a°6 (tc) 


J  s.gcw-o, 
an6(t) 


and,  by  (A3) , 


[  s  »g  cW  =  0 


B6nCt0t 


Further,  since  S  and  u  are  continuous  across  d  C.  ,  , 

’  ~o  ~o  o  t  t 

I  So'2od^=0’ 

B6nCt  t 
o 


a 
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+  _,1 


and,  since  s^eL  (C^)  and  u  is  bounded. 


(  s  *u  I  s„-u  dA'. 

j  ~o  ~  j  ~o  ~ 


B6nCt  t 
o 


V1 


Thus 


f  V  da  =  lim  J  V  da  =  J  (s-sQ)  •  (u-uq)  d#  -  J  sQ.u  d^. 


ft 


ft 


3ft 


'V 


where  we  have  used  the  fact  that,  by  (6.7)  and  (A^) ,  veL  (ft) . 

This  relation  also  holds  at  t  =  t  ,  since  (by  definition)  the 

integral  ov  ;r  is  zero.  Thus,  since 

o 

(at)t  I  =  J  { ^~s0)  *b  +  i-  io)  ]tU, 
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3ft 


0, 


t=t 


•  u 


dM 


this  relation  and  (6.6)  imply  (6.5)^. 

To  establish  the  equivalence  of  (6.5)^  and  (6.5) it  clearly 
suffices  to  prove  that 


lim 

W+oC 


J  Ste’V 


U(x,t)  -  u(X-t(t) ,t  ) 
f— - ^ - —IdA^ 


=  0 


t  t 
o 


Writing  x  =  zfc  +  r ,  we  see  that 

u(x,t)  =  u(zt+£,t), 
u(x-t(t),to)  =  u(zfc  +r,tQ)  . 


(6.8) 


30 


Thus  (cf.  the  paragraph  containing  (3.2))  the  quantity  in  {  ) 

in  (6.8)  is  bounded  in  absolute  value  by 

sup  lu7! , 
flx [0,T] 

which  is  finite  by  (A_)  .  Thus,  since  s-^eL*  (C,  ) ,  (6.8)  follows.  O 

z  T 
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